Robust Surface Doping of Bi 2 Se 3 by Rb Intercalation 
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Rubidium adsorption on the surface of the topological insulator Bi2Se3 is found to induce a strong 
downward band bending, leading to the appearance of a quantum-confined two dimensional electron 
gas states (2DEGs) in the conduction band. The 2DEGs shows a strong Rashba-type spin-orbit 
splitting, and it has previously been pointed out that this has relevance to nano-scale spintronics 
devices. The adsorption of Rb atoms, on the other hand, renders the surface very reactive and 
exposure to oxygen leads to a rapid degrading of the 2DEGs. We show that intercalating the 
Rb atoms, presumably into the van der Waals gaps in the quintuple layer structure of Bi 2 Se3, 
drastically reduces the surface reactivity while not affecting the promising electronic structure. The 
intercalation process is observed above room temperature and accelerated with increasing initial 
Rb coverage, an effect that is ascribed to the Coulomb interaction between the charged Rb ions. 
Coulomb repulsion is also thought to be responsible for a uniform distribution of Rb on the surface. 



INTRODUCTION 

The topologically guaranteed existence of metallic sur- 
face states on certain bulk insulators has recently at- 
tracted much interest in condensed matter physics [TJ [2] . 
While most research activities focus on the exotic prop- 
erties of the topological surface states, the study of the 
materials in question has also shown them to have highly 
promising properties for conventional spintronics applica- 
tions. Most importantly, strong band bending near the 
surface of the prototypical topological insulator Bi 2 Se3 
[H [4] , was shown to lead to quantized two-dimensional 
electron gas states (2DEGs) [5] with a very strong Rashba 
splitting [6] - something which could be used to realise 
nano-scale spin field effect transistors [7 j. The required 
band bending can be induced by the controlled adsorp- 
tion of a variety of species, such as CO |5], H2O [9], rare 
earth and alkali metals, [U[T0l[TT] but it is very sensitive 
to the amount of adsorbates, leading to an unstable elec- 
tronic situation. A stabilisation can be achieved for an 
adsorbate coverage that is sufficiently high to saturate 
the band bending process [SJ [10] but the particularly 
effective adsorption of alkali atoms can be expected to 
render the surface very reactive, even more than for the 
clean surfaces where this is already a problem p~2j [13] . 

Here we explore the possibility of intercalating ad- 
sorbed Rb atoms below the surface of Bi2Se3 and study 



what effect this has on the surface electronic structure 
and the surface stability. Compared with the situation of 
Rb adsorbates on the surface, we find that the intercala- 
tion of Rb atoms has almost no effect on the near-surface 
electronic structure as measured by angle-resolved pho- 
toemission spectroscopy (ARPES): the 2DEGs are al- 
ready formed when the Rb atoms are adsorbed on the 
surface and they are not severely changed by the interca- 
lation of Rb below the surface. Intercalating the highly 
reactive Rb atoms, on the other hand, protects them 
from interaction with the residual gas and the surface 
electronic structure is not even affected by a massive ex- 
posure to O2. Real space studies by scanning tunnelling 
microscopy (STM) show that one of the factors driving 
the intercalation process is the Coulomb repulsion be- 
tween the Rb atoms. This leads to a stronger tendency 
for intercalation at higher surface coverage and to a uni- 
form distribution of the Rb dopants. Such control of the 
dopant uniformity is advantageous if spin-transport in 
these systems is to be exploited. 



RESULTS AND DISCUSSION 

[l] shows the electronic structure of calcium-doped 
Bi2Se3 after cleaving and upon the adsorption of a sub- 
monolayer coverage of Rb at 190 K. Initially, the topo- 
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FIG. 1: (a) ARPES measurement of the near-surface elec- 
tronic structure of clean, calcium-doped Bi2Se3 through the 
Brillouin zone center and (b) the photoemission intensity at 
the Fermi level (high photoemission intensity is bright), (c) 
The corresponding electronic structure and (d) Fermi surface 
after deposition of sub- monolayer Rb. Measurements are per- 
formed using photon energies of 16 eV and 32 eV for (a),(b) 
and (c),(d), respectively. 



logical surface state with its characteristic Dirac cone is 
observed to be very close to the Fermi level, resulting in 
a small circular Fermi surface. After Rb deposition, two 
Rashba-split 2DEGs can be observed in the conduction 
band region, leading to a complex Fermi contour con- 
sisting of five concentric features, four from the 2DEGs 
and one from the topological surface state (two contours 
nearly coincide), consistent with previous observations 
[TO1 [TT] . The center contour is almost circular but a pro- 
gressive hexagonal warping is observed for the contours 
further away from the center, as expected from the bulk 
band structure [TJJ [15]. At the temperature of dosing 
(190 K) and measurement (70 K), it appears unlikely 
that the Rb atoms are intercalated below the surface but 
in the following we show that such an intercalation can 
be achieved by annealing the sample to higher tempera- 
tures, subsequent to low-temperature deposition. 

The temperature-induced intercalation of Rb atoms 
can be followed by monitoring the Bi 4f, Rb 3d and Se 3d 
core levels. The bottom spectra in[2]show the photoemis- 
sion intensity in the appropriate binding energy regions 
prior to Rb adsorption. The Bi 4f [16] and Se 3d core lev- 
els are each found to be well described by two doublets; a 
very small minority component which can be ascribed to 
surface adatoms, clusters or bound impurities and a ma- 
jority component assigned to atoms in the bulk of Bi 2 Se3. 
No surface core level shift for the outermost Se atoms is 
observed. The curve fits (black dotted lines) are super- 
posed on the data for both core levels (solid coloured 
lines). The intensity ratio and separation is fixed follow- 
ing literature values [171 [18]. In each case, the position 
of the largest component has been highlighted by a black 
vertical line. 

Once 0.23 monolayers (ML) of Rb are adsorbed at 
room temperature, a clear Rb 3d core level is also ob- 
servable. It can also be well described by two doublets 
(also shown in [2| , each with the expected separation and 
intensity ratio. The intensity of the low binding energy 
doublet (outlined in pink) is very weak. This component 
will subsequently be shown to arise from the intercalated 
Rb atoms and its presence immediately after adsorption 
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FIG. 2: (a) Bi 4f, Se 3d and Rb 3d core levels (coloured lines) 
and fits (black dotted lines) of a clean surface, following depo- 
sition of 0.23 ML Rb and incremental annealing for 5 minutes 
at the temperatures given. The Rb 3d core levels can be fitted 
by two doublets that are assigned to on-surface Rb (purple) 
and intercalated Rb (pink). Measurements are collected at 
room temperature using a photon energy of 353 eV and have 
been normalised such that the principal components are of 
equal intensity, (b) Relative intensity of the surface and inter- 
calated Rb components at selected stages of the experiment. 



suggests that a degree of intercalation occurs already at 
room temperature. It is also apparent that the adsorp- 
tion of Rb immediately leads to a shift of the Bi and Se 
core level lines by 0.4 eV to higher binding energy. This 
is the expected behaviour due to surface doping, con- 
sistent with a similar shift of the valence band features. 
Due to the electrostatic nature of the shift, both the Bi 
and Se core level lines change their position by the same 
amount. Except for this shift, the shape of these core 
levels is essentially unchanged. 

Annealing the sample to increasingly higher tempera- 
tures leads to the following changes in the spectra: the 
low binding energy doublet in the Rb 3d peak strongly 
increases at the expense of the high binding energy dou- 
blet (see |2jb)). We interpret this as being caused by 
the conversion of on-surface Rb atoms to Rb atoms in- 
tercalated below the surface, most likely in the van der 
Waals gaps between quintuple layers. After annealing 
to 417 K, the high binding energy doublet has almost 
vanished, suggesting that the majority of Rb atoms are 
now found below the surface. The relative intensity of 
the Rb 3d core level with respect to Bi 4f and Se 3d also 
decreases during the annealing. This decrease, as well as 
a higher sensitivity to on-surface Rb atoms observed in 
grazing emission measurements, is consistent with Rb in- 
tercalation to depths which are significant relative to the 
electron mean free path of ~ 6 A. A final change during 
annealing is a small reversal of the initial doping: the Bi 
4f and Se 3d peaks move slightly back to lower binding 
energies, as seen by following the solid lines that mark 
the position of the highest peak. 

Although the effect of Rb intercalation is very pro- 



3 




0.1 0.0 0.1 

k x (A" 1 ) 



FIG. 3: ARPES measurements of Ca-doped Bi 2 Se3 (a) after 
deposition of 0.23 ML Rb at 190 K and annealing for one 
minute to 350 K and (b), (c) and (d) following exposure to 
640 L, 1600 L and 6400 L of 2 , respectively, (e) A similar 
sample after a similar deposition of Rb at 190 K, but without 
subsequent annealing and (f) and (g), the same not-annealed 
sample following exposure to 110 L and 290 L of O2, respec- 
tively. Data are collected at 70 K using a photon energy of 
16 eV. In each case, a red dotted line indicates the approx- 
imate position of the first Rashba-split 2DEG, such that a 
comparison of the extent of doping is easier. 



nounced in the core levels, it hardly affects the electronic 
structure in the valence and conduction bands. This 
can be seen by comparing [3^a) with [3^e) where a sim- 
ilar amount of Rb has been dosed at 190 K, but the 
sample in[3ja) has been annealed (in this case to 350 K) 
whereas the sample in[3je) has not. In the annealed case, 
the electronic structure is slightly sharper, and the dop- 
ing strength is slightly less, but otherwise the two cases 
are remarkably similar. This is significant since the loca- 
tion of the Rb is entirely different in these two cases, yet 
the electronic structure and surface doping are essentially 
unaffected. 

While the annealing and intercalation hardly affects 
the electronic structure of the surface, [3] illustrates that 
the annealed surface is much more stable towards a reac- 
tion with oxygen: after 6400 L of O2 exposure, the dop- 
ing is marginally reduced to ^0.60 eV, but neither the 
2DEGs nor the topological surface state are significantly 
degraded. Conversely, a similar experiment in which Rb 
is dosed at 190 K without further annealing reveals that 
the surface is quickly degenerated by exposure to O2; al- 
ready after ^300 L, the conduction band quantum well 
states are barely discernible. 

The intercalation process and its origin can be fur- 



ther studied by STM.ga)-(d) show STM topographs of 
samples with two different Rb coverages before and after 
annealing at 400 K for 10 min. In all samples the major 
features are single, well separated Rb atoms, which ap- 
pear as circular protrusions with a height of ^150 pm. 
They are thus easily discernible from the trigonal features 
of Se vacancies or from sub-surface Ca dopants. [19]. In 
fact, the corrugation of the latter features is more than 
ten times smaller so that they can only be observed in 
images taken for a substantially lower Rb coverage. For 
the sample covered with 0.025 ML Rb (a,b), a few dimers 
are visible before annealing [see arrow in (a)] and after 
annealing some clusters are formed, most probably from 
residual gas adsorption [see arrow in (b)]. After anneal- 
ing, the coverage of on-surface Rb atoms is reduced by 
20% and the dimers have disappeared. In contrast, for 
the sample covered with 0.12 ML Rb (|4]c,d), which is 
closer to the coverage used in the ARPES and core level 
photoemission measurements, there is a much stronger 
effect of the annealing on the distribution and coverage 
of on-surface Rb atoms. |4jc) shows the STM topograph 
of the non-annealed sample after Rb adsorption. We find 
a distribution of well separated Rb atoms with an in- 
crease in the apparent height in areas of locally higher 
coverage, probably due to charge accumulation effects 
[20] |2T] . After annealing, the coverage of on-surface Rb 
atoms is strongly reduced to 0.05 ML Qd)) indicating 
that 60% of Rb has diffused into the bulk. Moreover, the 
distribution of the Rb remaining adsorbed on the sur- 
face is much more homogeneous after annealing. The 
enhanced tendency for intercalation at high coverage, 
as well as the trend for a more uniform distribution of 
the Rb atoms after annealing, can both be ascribed to 
the strong Coulomb interaction between the Rb atoms 
that are highly ionic due to a charge transfer to the bulk 

P3JI22]. 

Our results show that the annealing of the Bi2Se3 sam- 
ple after adsorption of Rb leads to a decrease of Rb cov- 
erage on the surface. While we do not directly show that 
the disappearing Rb atoms are intercalated into the van 
der Waals gaps between the quintuple layers making up 
the Bi2Se3 structure, this appears very likely. In fact, 
such intercalation processes have been known for a long 
time [23] and intercalation can be achieved during bulk 
crystal growth [24] or using electrochemical methods [25] . 

Disappearing on-surface Rb atoms after annealing 
could be explained by thermal desorption instead of in- 
tercalation, as has been proposed in the case of potassium 
doping. [10] However, this alternative mechanism cannot 
be reconciled with the other experimental findings in the 
present case. Most importantly, a thermal desorption of 
the Rb atoms could be expected to reverse the observed 
band bending, such that the initial electronic structure is 
recovered, something that clearly is not the case. Desorp- 
tion would also lead to a mere decrease of the Rb 3d core 
level intensity, not to the appearance of a new component 
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FIG. 4: (a,b) STM topograph of 0.025 ML Rb on Ca-doped 
Bi 2 Se3 before (a) and after (b) annealing at 400 K for 10 min. 
After annealing the on-surface Rb coverage is reduced by 20% 
(b). Arrows: see text. (c,d) Corresponding data for an initial 
coverage of 0.12 ML Rb before (c) and after (d) annealing. 
After annealing the on-surface Rb coverage is reduced by 60% 
(d). The insets (10 nmxlO nm) show magnified views of 
the samples. (Tunneling parameters: V =1 V, I =15 pA, 
measurement temperature: 1.2 K.) 



in the spectrum. Finally, the desorption of alkali atoms 
from metal [26] or semiconductor [27] surfaces typically 
takes place at much higher temperatures. Since thermal 
desorption appears not to play a role here it is somewhat 
surprising that it should for the chemically similar K. 
Further direct structural investigations, for example by 
STM, should be able to clarify this point. An interesting 
similarity between thermal desorption and the intercala- 
tion observed here is the fact that higher alkali coverages 
favour both processes. In the case of desorption, a higher 
coverage lowers the desorption temperature, here it low- 
ers the temperature required to achieve intercalation. In 
thermal desorption, this effect can be described by a sim- 
ple model that takes the Coulomb interaction between 
the adsorbates into account [26] . 

So far, the intercalation of dopants from surface adsor- 
bates has only been reported for the case of Ag on Bi 2 Se3 
[28] , where this is observed already at room temperature. 
Also in that case, the Ag adsorption leads to a shift of 
the entire electronic structure to higher binding energies 
and to the appearance of additional states in the con- 
duction band region. This, however, is not interpreted 
as band bending but ascribed to an increased spacing of 
the van der Waals gap, caused by the intercalation of 
Ag. Indeed, first principles calculations can bring about 
features similar to the experimental ones, but only if an 
increase of the van der Waals gap by 36% is assumed. A 



similar increase of the van der Waals gap had also been 
evoked previously to explain not only the Rashba-split 
2DEGs in the conduction band but also the M-shaped 
features in the valence band [29 that emerge after the 
adsorption of Rb (see|3|. 

Band bending, on the other hand, can not only explain 
the formation of the Rashba-split 2DEG in the conduc- 
tion band region but also the M-shaped states in the 
valence band. These can also be described as quantum 
well states caused by an overall downward bending of the 
bands. In the case of valence band states, this is unusual 
because only upward band bending would be expected 
to lead to quantized hole states. In the special case of 
Bi2Se3, however, it can also quantize the valence band 
because the upper valence band is rather narrow and 
placed over a projected bulk band gap [8 . In this pic- 
ture, a full quantization of the valence band states would 
require a band bending larger than the total bandwidth 
300 meV) but the lowest lying states could already be 
quantized for a smaller band bending. 

Evidently, Rb adsorption has a very pronounced effect 
on the near-surface electronic structure but the formation 
of 2DEGs in the conduction band and of the M-shaped 
states in valence band region already takes place when Rb 
atoms are adsorbed on the surface. The further interca- 
lation brings about only very minor additional changes. 
On the other hand, an increased van der Waals interlayer 
spacing upon intercalation is very likely but it is probably 
quite small. In fact, bulk x-ray diffraction measurements 
have found such an increase for Cu-intercalated Bi 2 Se3 
to be smaller than 3% [24] . 

Another important result is that the intercalation of 
Rb not only increases the chemical stability, for exam- 
ple towards oxidation. The saturated strong band bend- 
ing also results in a very stable situation towards further 
changes in the electronic structure by the adsorption of 
other species that otherwise increase the band bending. 
The electronic structure of clean Bi2Se3 is very sensitive 
to the adsorption of minute quantities of residual gas in 
ultra-high vacuum (UHV) [30 because this causes a band 
bending. The adsorption of Rb improves this situation 
because the band bending can be saturated, similar to 
the situation recently reported for K adsorption at low 
temperature [TO] . Exposing the non-intercalated system 
to UHV for 24 hours leads to a reduction of the band 
bending by ~ 100 meV. The full chemical protection is 
achieved by the intercalation of Rb that results in a sur- 
face that remains stable in UHV overs several days. 



CONCLUSION 

We have demonstrated that Rb atoms adsorbed on 
Bi2Se3 can be intercalated below the surface by anneal- 
ing above room temperature. This results in a very stable 
situation with respect to further band bending and chem- 
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ical reactivity. Intercalation of adsorbates is a promising 
way to influence the electronic structures of bulk topolog- 
ical insulators where it is most relevant: near the surface. 
Controlled intercalation processes could be an important 
step in the fabrication of topological insulator-based de- 
vices before further surface passivation is carried out [13] . 
In particular, the intercalation of species giving rise to an 
upward band bending could lead to a space-charge layer 
that is depleted of electrons, even if the bulk material is 
strongly n-doped, as often the case for these materials. 

EXPERIMENTAL DETAILS 

Bulk Bi2Se3 crystals were calcium doped to position 
the bulk Fermi level within the bulk band gap. These 
samples were synthesized by melting a mixture of Bi (5N 
purity) , Se (5N purity) and Ca (99.5% purity) with ratio 
Bi : Se : Ca = 1.996 : 3 : 0.004 at 860° C for 24 hours in 
an evacuated quartz ampoule. Once cooled down from 
860° C to 750° C at a rate of 50°C/h and then from 750° C 
to 600° C at a rate of 2 °C/h, the sample was annealed 
at 600° C for 7 days. 

ARPES and core level photoemission measurements 
were performed at beamlines SGM III [31 and SGM I, re- 
spectively, at the synchrotron radiation source ASTRID. 
Clean Bi2Se3 surfaces were obtained by cleaving bulk 
crystals in UHV, and the surface state Dirac point was 
measured to be less than 50 meV below the Fermi level. 
Exposure to Rb and O2 as well as annealing were per- 
formed in situ. The surface coverage of Rb was deter- 
mined by the relative intensity of Rb and substrate core 
level lines. It is given in monolayers (ML), i.e. relative 
to the number of Se atoms in the first layer. The sur- 
face was exposed to molecular oxygen and the exposure 
is given in Langmuir (L). Rb was typically dosed within 
20 minutes of cleaving, such that surface doping due to 
ageing can be considered insignificant. 

We have also carried out ARPES experiments for Rb 
adsorption on pristine Bi 2 Se3 samples to test the the role 
of the Ca doping. As expected due to the small quantity 
of Ca, we have found no indication of different behav- 
ior. Annealing the Rb-dosed samples, one can achieve a 
smaller but equally stable band bending and a surface 
that is stable in UHV over several days. 

The STM experiments were performed in an UHV 
low-temperature scanning tunneling microscope facility 
cooled by a Joule-Thomson cryostat [32]. The bulk 
Bi2Se3 samples were cleaved in UHV at room temper- 
ature and immediately transferred into the cryogenic mi- 
croscope. During cool down, Rb was deposited onto the 
substrates at T j 200 K. For annealing, the samples were 
transferred in UHV from the STM to a Boron Nitride 
heater with a thermocouple. STM topographs were ac- 
quired at 1.2 K in the constant current mode at a current 
/ with a bias voltage V applied to the sample. The cov- 



erage of on-surface Rb was determined by counting the 
number of Rb atoms visible in STM topographs of a cer- 
tain area and dividing by the number of surface Se unit 
cells within the same area. We have also performed an- 
nealing experiments for the bare calcium-doped samples. 
For annealing temperatures below 423 K, we did not no- 
tice segregation of Se or Ca to the surface. 

The use of three separate experimental systems means 
that there can be a systematic uncertainty between the 
temperatures quoted, which we estimate could deviate 
by up to 20 K from one system to another. It also means 
that experimental parameters such as Rb dose, and an- 
nealing and measuring temperatures, are not identical. 
Since similar behaviour is observed over an extended pa- 
rameter range, this is not thought to be significant. Also, 
the Rb coverages are determined with a rather different 
accuracy in the STM and the photoemission setups. The 
coverage determination by STM is highly accurate but 
the uncertainty in the photoemission measurements can 
be larger, and may be up to 20% for the quoted values. 
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